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scattering of 92.5 MeV incident electrons at scattering 
angles of 75°, 90° , 105° and 120°. In the excitation energy 
range of 6.1 to 38 MeV, nine transitions were observed. 

The previously reported El has been verified, but the broad 
E2 has been separated into two distinct resonances. Reduced 
transition probabilities and multipolarity assignments have 
been made. Resonances occurred at excitation energies of 
6 . 69 (E2) , 8.09 (E2) , 10.01 (E2) , 11.21 (E2) , 12.46 (E3) , 
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I. 



INTRODUCTION 



Inelastic electron scattering has become the major tool 

of investigation of giant resonance phenomena. Previous 

experiments with the Naval Postgraduate School Linear Accel- 

208 197 

erator have included Pb, Au and Ho [Refs. (1) and 

(2)]. In order to make a more systematic survey of the 

nuclear table, it was decided to investigate medium-light 
. 89 

nuclei. Y was chosen because it has a closed neutron 

shell (N=50) . Previous experience with closed shell nuclei 

[Refs. (3) and (4)] have shown that giant resonances in these 

nuclei are relatively small and can therefore be disentangled 

. 90 

from each other. Furthermore, the neighbor nucleus Zr 
is the only nucleus where a total E2 width (4.8MeV) greater 
than the total El width (4.0MeV) has been reported [Ref. 

(4) ] . 

With a fixed energy incident electron beam of 92.5 MeV, 
data were collected at scattering angles of 75°, 90°, 105° 
and 120°. For each scattering angle, the number of scattered 
electrons was measured in the range 96 MeV to 50 MeV in order 
to include the elastic peak and excitation resonances from 
6 to 38 MeV. Experimental inelastic cross sections were 
extracted for nine states and multipolarities have been 
assigned . 
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II. THEORY 



A. ELECTRON SCATTERING 

1. Elastic Electron Scattering 

The principle of electron scattering from nuclei 
can be found in the derivation of scattering theory cross- 
sections by Rutherford and Mott. The Rutherford differential 
cross section for elastic scattering of charged point parti- 
cles from a fixed point charge Ze is written [Ref. (5)]: 



(— ) 

v dft' RUTH 



16 E i 2 sin 4 |- 



(II-l) 



• where 



9 = scattering angle 
E^ = total energy of incident particle 



This equation does not account for relativistic and spin 
effects. Mott, following Dirac's formalism for the relati- 
vistic electron, took the latter effects into account. In 
considering scattering from point charge centers, Mott 
derived the following formula for extremely relativistic 
Dirac particles [Ref. (6)]: 



(— ) 

v dJT MOTT 



Z e‘ 



2E 



- 2 9 

cos 7 

. 4 9 

Sin 



(II-2) 
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This formula was followed by experimental confirmation 
using electron energies on the order of 1 MeV [Ref. (7)]. 

In more recent years, electron scattering has been trans- 
formed into the rapidly evolving science of nuclear structure 
study. Using electron beams of sufficiently high energy 
(E 50 MeV) so that the electron's De Broglie wave length 
(X = h/k) becomes equal to, or smaller than, the spatial 
extension of the nucleus ( - 1FERMI = 10 ^m) , researchers 
began to probe the structure of the nucleus [Ref. (7) ] . 

The angular distribution of electrons scattered is influenced 
by the extension and shape of the nuclear charge distribution. 
Therefore, experimental results which deviate from Mott 
cross section predictions are interpreted as arising from 
the finite extent of the nuclear charge density. To account 
for the finite size of the nucleus, a form factor representing 
the charge distribution of the target nucleus multiplies the 
Mott cross section [Ref. (8)], giving 

( dfi } = (d^MOTT l F(c J ) l ' (H-3) 



where 



q = momentum transfer, and 
F (q) = / p (r) e iq ’ r d 3 r. 

A minor correction must be applied to Equation II-3 to account 
for recoil of the target nucleus [Ref. (8)]. 
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Application of these principles to elastic electron 
scattering has yielded a great deal of information about 
both nucleon and nuclear ground state charge structure. 
However, this knowledge represents only a small portion 
of possible nuclear information since it bears only on the 
static properties of the ground state of the nucleus. The 
prospect of exciting the nucleus to higher energy levels 
in order to observe nuclear dynamics led to the study of 
inelastic electron scattering. 

2 . Inelastic Electron Scattering 

Inelastic electron scattering is one of the most 
powerful tools available for nuclear structure studies. 

There are two basic reasons for this. First, the basic 
interaction is the well known electromagnetic interaction 
of the electron with the nuclear charge and current. 

Second, and more important, is the possibility of separately 
varying energy and momentum transfer (q) , allowing mapping 
of the inelastic form factors of nuclear levels as a func- 
tion of the momentum transfer [Ref. (7) ] . 

Elastic scattering deals with ground state proper- 
ties only. Inelastic scattering involves the excitation from 
the ground state to various excited states. When a rela- 
tivistic, monoenergetic electron beam is incident on a thin 
(t << radiation length) target, a small fraction of the 
electrons will undergo collisions with nuclei. Use of a 
magnetic spectrometer capable of rotation about the scattering 
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axis enables the determination of the energy distribution 
of the scattered electrons at various angles. Most electrons 
are scattered elastically and appear in the spectrum as a 
sharp peak at an energy lower than the beam energy due to 
recoil, and with a width determined by overall experimental 
resolution. Below the elastic peak there appears in the 
spectrum of scattered electrons a continuous distribution 
called the radiation tail, which is due to bremsstrahlung . 
Superimposed on the radiation tail are peaks associated with 
nuclear excitation. All of the scattered electrons below 
the elastic peak have lost energy due to one or more energy 
exchange mechanisms. The experimentalist is only interested 
in those events involving the exchange of one photon from 
an electron to a nucleus, resulting in excitation of the 
nucleus. Therefore, all other energy exchange mechanisms 
must be understood well enough to accurately correct the 
inelastic spectra for proper observation of the superimposed 
nuclear resonances. 

These corrections are most important for the elastic 
scattering, but they apply similarly to inelastic excitation. 
They fall into two categories: line shape corrections and 

radiative tail corrections . Line shape corrections arise 
since some of the scattered electrons are degraded in 
energy prior to being counted by the emission of soft pho- 
tons, emission of a hard photon (bremsstrahlung), multiple 
collisions with nuclei, collisions with atomic electrons and 
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straggling due to ionization effects. These all lead to 
multiplicative factors, called radiative corrections, applied 
to the area under the elastic peak, which is proportional to 
the elastic cross section. The radiative tail is the energy 
distribution of the electrons which, by the processes just 
described, have been scattered out of the energy region of 
the elastic peak and form a continuous spectrum under the 
inelastic levels [Ref. (7)]. Naturally, each inelastic 
level has its own radiation tail which, for practical pur- 
poses, is mostly neglected in the evaluation. 

Once the radiation tail and other experimental 
background have been subtracted, the remaining cross section 
consists of the inelastic cross sections corresponding to 
nuclear excitation levels. The analysis of how these cross 
sections depend on the scattering angle and therefore on 
momentum transfer q, gives the experimentalist a tool for 
studying nuclear structure and dynamics. 

Since various physical parameters of the experimental 
arrangement are not known with sufficient accuracy, it has 
become customary in inelastic scattering to determine the 
ratio of inelastic to elastic cross section by measuring the 
inelastic peak area (A^ n ) relative to that of the elastic 
peak area (A^) , instead of doing absolute measurements. 

For targets of pure isotopes, this area ratio is proportional 
to the ratio of corresponding cross sections. 
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(II-4) 




where K is a correction factor which accounts for: 1. dif- 

ferences in the radiative and ionization corrections of the 
elastic and inelastic peaks; 2. double or multiple scattering 
(more than one photon exchange); 3. the result of averaging 
the cross section over the spectrometer acceptance solid 
angle; and 4. any apparatus effects. Thus, the accuracy 
to which inelastic cross sections can be determined is 
limited by the experimental uncertainties of the peak area 
values and the precision to which the elastic cross section 
is known [Ref. (9)]. 



F(q) for elastic scattering from the theoretical and experi- 
mental cross sections. An analogous (but more complex) 
equation can be developed for inelastic scattering. The 
form factors in these equations contain information on the 
nuclear ground state for elastic scattering and on a given 
nuclear excitation level for inelastic scattering. To see 
this quantity more clearly, division of both sides of 
equation (II-4) by the Mott cross section yields 



Equation (II-3) provides the nuclear form factor 



(— ) • 
'dJr in 



c ( dft } el 



df2 MOTT 



'dfl'MOTT 



where: 
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A. 

C = K (~i) 
el 



The use of such a relative cross section eliminates most 
of the kinematic factors, makes nuclear contributions more 
evident, and allows for more convenient presentation of 
the data in analogy to that of elastic scattering [Ref. 

( 10 ) ] . 

B. DISTORTED WAVE BORN APPROXIMATION 
1. Cross Section Calculation 

Actual computation of the cross section has become 
a more sophisticated procedure in recent years. Originally, 
the incoming and outgoing electron was treated in the Dirac 
formalism as a plane wave. In the Born approximation (or 
Plane Wave Born Approximation) , the differential cross 
section can be written as a sum over the cross sections for 
electric (E) and Magnetic (M) multipole transitions [Ref. 

( 9 ) ] : 



/da n 

( an> 



PWBA 



= Z 
X 



(£2.) 

W 



+ E 



EX 



(d£) 

W 



MX 



(II-6) 



where 



( a§>Ex = a 2 a x q 2 V 2 w + 1 ) ‘ lB(cX ' q ' I o* I x >v L (9) 



+ B(EX,q,I 0 +I x )V T (e) ] R 



-1 



and 
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( 3 §> . “ C‘ 2 a x q 2 X k 0 ' 2 B(MX,q,I 0 *I x )V T (e)R - 1 . 

In the above 

a x = 4 itX~ 1 (X+1) [ (2X+1) II]” 2 

k o = V Hc 

R = 1 + hc(k o /mc 2 ) (1 - cos 9) 

X = transition multipolarity 
a = 1/137 = fine structure constant 
E q = primary electron energy 
9 = scattering angle 

m = nuclear mass. 

However, steadily improving accuracy in electron 
scattering experiments has required more accurate analysis 
of data to obtain meaningful results. By modeling the 
incoming and outgoing electron as waves distorted by the 
Coulomb potential of the nucleus, this has been accomplished. 

An early publication of such an analysis dealing with elastic 
scattering was presented by Yennie and Ravenhall in 1954 
[Ref. (11) ] . Ravenhall used the model of the Dirac equation 
for an extremely relativistic electron in the elctrostatic 
potential of a static spherically symmetric charge distribution. 
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and found a new cross section obtained by a numerical 

calculation of the phase shift of each partial wave. 

Electron scattering is concerned with high energy 

electrons (E^ > 30 MeV) where the electron rest energy can 

2 

be neglected (me is effectively zero). Therefore, the 
Dirac equation can be written as 

(a -pc + V - E^S = 0 , (II-7) 

for which the scattering states have the asymptotic form 

$ ~ (J) e lkz + r -1 f ( 0 , <f>) ( q ±( J e lkr (II-8) 

tan 2 e ^ 

which is the same for both spin orientations when V is a 
spherically symmetric potential. This can be decomposed 
into partial waves and summed: 



$ 



Z 





(II-9) 



where: 



J 



2 




j ( j +1) h 2 




and 



J 



z 




m h 

]m 
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with 



"t — * 1 7“ 

J=rxp + jah . 



By applying scattering theory at large distances 
from the scattering center and modifying the results to 
account for the long range effects of the Coulomb potential, 
the scattered wave becomes 



*SCATT - X ' 1 k f < 9 > < 9 i*> e i(x+ vln2 X ) (II . 10) 

tan 2 e * 



where 



y = Z 



e 2 /hc 



and 



f ( 0 ) = 



1 

2ik 



I e 



2in j 



(j + i) (P j . % + P j+h ) 



for 



e t* o. 

However, the phase shifts ( n j ) do not approach a 
limit as j increases and therefore the summation of f ( 0 ) 
is difficult and in practice is done using computer codes 
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[Ref. (11) ] . Later, Rawitscher and Fischer modified the 
calculation to apply to cases where the energy is not in the 
extreme relativistic range [Ref. (12) ] . 



Application of these principles to inelastic 



scattering is more involved and proceeded more slowly. The 
first successful computer codes for solving equations in the 
inelastic case was called GBROW and was described in a report 
by J. F. Ziegler [Ref. (13)]. Five conditions were assumed 
in GBROW: 

1. Exchange of a single photon. 

2. The transition is of a pure electric multipole character. 

3. The nuclear ground state is spherically symmetric. 

4. The excited state charge distribution is not signifi- 
cantly distorted from the ground state. 

5. Nuclear recoil is negligible. 

In this program the inelastic cross section is derived in 
the form of Fermi's Golden Rule: 



(— ) • 
v d£r mel 



(kinematic terms) £ 




i 



( 11 - 11 ) 



where 



H = static interaction hamiltonian 




o 



H^ nt = dynamic interaction hamiltonian 



and 
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/ 



fi over nuclear volume f 



ip T (N) 



(electron field) 



ip (N) 



(N) 

The integration is over the nuclear volume and \p ^ and 
(N) 

\p ^ are the final and initial nuclear wave functions. 

Here the "electron field" is an expansion of the electron's 
electromagnetic field in partial waves. 

Proceeding to separate Coulomb and current portions of 
the Dirac particle interaction, orienting the electron Z 
coordinate axis parallel to the incident electron direction, 
and using orthogonality relations between various spherical 
harmonic terms, the resultant cross section is 



(— ) . 
v dJT inel 

DWBA 



V 2 E i E f P f ( 2 I f + 1) 

2(2„ 2 )K 3 P i T2X-TTTTlL7Tr <H i nt> fi 



(II-12) 



where 



V = normalization volume, 

L = transition multipolarity, 

and 



<H. .> is defined in [Ref. (13)]. 

mt fi 

If the cross section is written as in Equation (II-3) , the 
form factor squared has the form (in Plane Wave Born 
Approximation) 
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F l 2 = ( 1 p c ,lI 2 + + tan2 l )[ |f E , L l 2 + l F „, L l 2 n (h-13) 



2 . 



where | F c L | is due to the Coulomb interaction only and 
is written 



F 



C,L 



2 



1 4 it _2L 

“T 9 q 

Z z [(2L+1)!!] 



B (cL,q, J i 



J f ) 



(11-14) 



where 



Z = atomic number 
q = momentum transfer 

= initial/final angular momentum of nucleus. 



The transverse electric and magnetic terms are 



E ,L 
M 



4tt ,L+1. , q . 2 _ ,E t t t . 

2 L *• ( 2L+1) ! ! 1 B ( M L ' q ' 

z 



(11-15) 



The coefficients B(X) in the above are the reduced nuclear 
transition probabilities. They now contain the desired 
information on a given level and are of the form 



B(XL f q,J i ^J f ) = 2J 1 +1 |<J f | |M(XL,q) | |J i >| 2 , (11-16) 

i 

A 

where M represents the transition operator and will, in 
practice, contain a function representing the nuclear model 
used. In the distorted wave treatment, this function is the 
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nuclear model-dependent transition charge density operator 

P tr (r). Determination of spacial charge density functions 

(r) , p 2 (r) , p 3 (r) from this operator will then give an 

evaluation of the reduced transition matrix elements and 

enable the cross section calculation to be completed [Ref. 

(13)]. The model used in GBROW is the hydrodynamic model 

developed by Tassie [Ref. (14) ] . The B- value is then 

2 

defined and by B (EL) = /p tr (r) r dx for electrical 
longitudinal transitions. 

2 . Tassie Model 

It is customary to approach the actual calculation 
of the DWBA cross-section by first choosing an appropriate 
nuclear model and its associated transition charge density 
P tr (r) and transition currents j N » The model is considered 
satisfactory if it furnishes computed cross section values 
that agree with the data. It has been found that for collec- 
tive transitions (many particle excitation) , reasonable 
agreement is obtained with the use of an oscillating liquid 
drop model described by Tassie [Ref. (14) ] . This model is 
based on hydrodynamic vibrations of the nucleus. For collec- 
tive excitations, transverse contributions to the cross- 
section are small below angles of about 150°, and can there- 
fore be ignored without much loss of accuracy [Ref. (10) and 
(15)]. As presented by J. F. Ziegler in 1967 [Ref. (13)], 
the Tassie model transition charge density for a multipolarity 
L is [Ref. (14) ] 
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(H-17) 



P L = r N 



L-l ^o 
3r N 



This form of p is then entered into the DWBA calculations 
of the cross-section by the computer code GBROW. 

C . GIANT RESONANCES 

Prior to the development of inelastic electron scattering 
techniques, electromagnetic excitation of nuclear levels 
was studied mainly through photonuclear methods. A dominant 
feature of these spectra was known for years as "The Giant 
Resonance" or "Giant Dipole Resonance." This structure con- 
sistently appeared in spectra throughout the periodic table 
for excitation energies between 10 MeV and 25 MeV. Its 
absorption cross section exhausts approximately all of the 
classical Thomas-Reiche-Kuhn (TRK) dipole sum rule, which 
was originally derived for atomic excitations. It was 
therefore regarded as a strongly collective excitation in 
which a considerable number of nucleons participated. 

Through the use of monocromatic beams, produced by positron 
annihilation, very accurate results for excitation energy, 
width and line shape of the GDR were obtained [Ref. (16)]. 

During this period, electron scattering experiments 
were concerned mainly with investigations of the nuclear 
ground state, low level excitations and back scattering 
(180°) to isolate magnetic transitions, although some GDR 
studies were done in light nuclei. With the advent of 
better facilities, the investigation of giant resonances 
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in the nuclear continuum became possible. The first measure- 
ments covering nuclear excitations up to 28 MeV showed not 

140 

only the GDR(El) resonance in Ce but also an E2 giant 
quadrupole resonance (GQR) below the GDR at 12.0 MeV which 
accounted for 65% of the E2 isoscalar sum rule and a giant 
magnetic dipole state (Ml) at 8.7 MeV [Ref.s (3), (17)]. 

One of the characteristics of giant resonances is the 
relatively smooth dependence of their excitation energies 
(E ) with the mass number, A. The use of a simple hydro- 

A 

-1/3 

dynamic nuclear model yields predictions of E (GDR) = 80A 

A 

for the GDR, close to the experimentally observed value in 
heavy nuclei. The GQR (E2) had been predicted earlier by 
Bohr and Mottelson [Ref. (18) ] as an isoscalar collective 
E2 mode. Their predictions were 

E (E2) = — (2hw ) = 58A (for the isoscalar E2) 

x /2 ° 



and 



E (E2) = 135A (for the isovector E2) . 

X 

140 

The isoscalar prediction of 11.2 MeV m Ce compares 
favorably with the Darmstadt results of 12 MeV [Ref. (19)]. 

The discovery of the "new" giant resonances (GQR) in 
inelastic electron spectra results from the fact that the 
momentum transfer, q, can be varied independently of 
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in electron scattering. The y~ 



excitation energy, E , 

X 

absorption cross section, for a given energy, can be 
measured only at q = E/c, which is called the photon point. 
However, in electron scattering, q can be varied over a wide 
range above the electron energy line by varying either the 
primary energy or the scattering angle, or both [Ref. (7) ] . 
The governing equation for q(E,6) is 



q 



2 



4 E i E f 



. 2 9 
sm j , 



(11-18) 



where 



E^yf = electron initial/final energy. 

By thus varying q, giant resonances with different 
multipolarities having different q dependence can be distin- 
guished. Each of these resonances is associated with a 
particular value of angular momentum, L, transferred to the 
nucleus during the excitation. The GDR is an El mode 
corresponding to an L = 1 or dipole excitation. Experiments 
in many nuclei have produced evidence of EO, E2, E3 and Ml 
resonances within the nuclear continuum. 

Q Q 

D. Y EXPERIMENTAL SURVEY 

For comparison to the present measurements, a survey 

89 

was made of previous experiments with Y, for which there 
have been fewer experiments than with other nuclei in this 
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8 8 

mass range. Initial efforts were directed towards an Sr 

core with an extra proton [Refs. (20, (21) and (22)] but 

the weak coupling model could not adequately describe the 

excited states. Continued effort was almost strictly in 

the 0 to 4 MeV range and is tabulated by date in Table I. 

In the giant resonance region, a radiative proton cap- 
8 9 90 

ture reaction Y(p,y o ) Zr determined the existence of 

90 

giant El excitations in Zr at E = 16.5 MeV with a width 

of r = 4 MeV [Ref. (23)]. Using (e,e'), Fukuda and Torizuka 

9 0 

measured El at E = 16.65 MeV and T = 4.0 MeV in Zr. 

X 

The latter experiment also reported an E2 at E =14.0 MeV 

X 

and T = 4.8 MeV and at E = 27MeV. All fits were made 

X 

using Lorentz line shapes. Here, incident electrons of 150 
to 250 MeV were scattered at angles of 35° and 45° speci- 
fically so that the longitudinal terms would dominate the 

excitation function. Youngblood [Ref. (24) ] reported an 

90 

E2 giant resonance found in Zr using a particle scattering 
at E = 14.5 ± .3 with width V = 4.0 ± .2. This resonance 

A 

took up 56% ± 17% of the E2 sum rule. 

In an as yet unpublished paper [Ref. (25)], Bertrand 
has reported an isoscalar giant quadrupole resonance E2 

O Q 

in Y at E = 13.8 ± .2 MeV of width T = 3.32 MeV. These 

X 

results are from an (p,p') experiment and contain 24% ± 5% 
of the energy weighted sum rule. The measurements were taken 
at Oak Ridge National Laboratory . 
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A compilation of photonuclear experiments by Berman 
[Ref. (16)] gives data on the giant dipole resonance El 
at E = 16.77 MeV and r=4.1 MeV. This data compares favorably 

A 

with the present work. 
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TABLE I. PREVIOUSLY REPORTEC TRANSITION 



ENERGY 

(MEV) 


WICTF 

(MEV) 


TYFE 

EXPERIMENT 


REFER Ef*C E 


1.5X) 


RES 


( N » N ' Y) 


20 






(OC » ) 


26 






<P,P' ) 


21 






(E, E' ) 


22 


1.74 


RES 


(N,N» 


20 






(*, <*') 


26 






(P,P'J 


21 






( E , E • ) 


22 


2.21 


RES 


<n » n ' y ) 


20 






(P»P' ) 


21 






(E, E* 1 


22 


2.52 


RES 


( n » n ' y ) 


20 






( p , p • > 


21 






( E » E ' ) 


22 


2.86 


RES 


( N » N ' y ) 


20 






(E, E' ) 


22 




DOUBLET 


( P , P • ) 


27 


3.1 


RES 


( n t n * y ) 


20 






(E, E* ) 


22 




TRIPLET 


(Pt P' ) 


27 


2.72 


RES 


( E » E ' ) 


22 


4.0 


RES 


< E , E • 1 


22 


4.46 


RES 


( E t E ' ) 


22 


12. £+.2 


2.2 


(P, P* ) 


25 


16. 7S 


3.95 


(Y, Y’ ) 


28 


16 .74 


4.25 


( y » y‘) 


25 




NOTE: RES = RESOLUTION 
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III. DATA ACQUISITION 



A. EXPERIMENTAL PROCEDURES 

The physical layout and operational procedures of the 
NPS LINAC have been well documented [Refs. (30) and (31)]. 

A major improvement has occurred in the meantime with the 
installation of ion vacuum pumps to replace oil diffusion 
pumps in much of the vacuum system. Extensive acid cleaning 
of the wave guides to remove oil deposits was performed 
before the completion of the vacuum pump replacement. This 
change resulted in more reliable operation of the LINAC. 

Several runs were made prior to the actual data acquisi- 
tion in an effort to check and reduce background in the 
counting system. Of the four angles measured, the 105° 
measurement had the least background, because of extensive 
extra shielding added to the beam pipe just following the 
deflection magnets and prior to the target chamber. 
Fluctuations in the data were also reduced by improvements 
to the counting system beam current monitor. This device 
is designed to shut off the counting system whenever the 

analyzed beam current falls below a preset value. 

12 

Two C target calibration runs were made to calibrate 
the magnetic spectrometer by checking the measured excitation 
energy of the 15.1 MeV state. 
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B. DATA ACCUMULATION 



A total of seven experimental runs were made using 

92.5 MeV incident electrons scattering at 75°, 90°, 105° 

and 120° deflection angles. The scattered electron energies 

were measured from 96 to 50 MeV in 0.1 MeV steps. Targets 

2 2 

with a width of 180 mg/cm and 280 mg/cm were used. The 
procedures for all runs were identical to those described 
in Ref. [(31)]. 

C. DATA REDUCTION 

The single counter spectra for the ten counters of the 
counter ladder were collected by means of a teletype printer 
and tape punch unit. The tape was compiled and read onto 
magnetic tape by use of the NPS Digital Equipment Corporation 
PDP11-50 Duplex System. This tape was then read onto the 
NPS Computer Center IBM 360/67 time-sharing system for text 
editing and data reduction purposes. The final spectrum 
was then stored on the mass storage (Data Cell) of the OS 
System for further evaluation. Determination of the 
resonances which make up the inelastic spectrum was accom- 
plished with a line shape fit program. The elastic line 
shape was assumed to be made up of two Gaussian curves both 
with the same height but different widths, which were joined 
to a short radiation tail of assumed hyperbolic shape. The 
elastic cross section, thus calculated, was then used to 
calculate the radiation tail and inelastic cross section 
(see Section II. A. 2., eqs. II-4, 11-12). 
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As a result of the comparison of Breit-Wigner forms to 
Lorentz forms in Ref. [(32)], and the better fits experienced 
with the former, Breit-Wigner line shapes were used for the 
resonance fits. The photonuclear work described in Section 
II. D. gives the parameters for the El resonance, E =16.77 

X 

and T = 4.1 MeV, which were used as a starting point. The 
B-value of the El resonance can be calculated from the 
integrated cross section by using the formula [Ref. (15)] 

/ a dE = 7r 2 Eca — — ■ — ■ 9 k 2L_1 B (L,k) . (III-l) 

Y [ ( 2L+1) ! ! j 

The peak height of the El resonance in the different spectra 
was then chosen to reproduce the calculated B-value of 
19.7 fm 2 . 

The following criteria to determine a reasonable fit were 
used: 

(1) The data and calculated spectrum should coincide 
visually . 

2 

(2) The x P er degree of freedom should be less than 
one. The errors are not strictly statistical because the 
detector momentum interval is larger than the momentum 
increment of the spectrometer field and hence correlations 
exist between energy bins. 

(3) All observed resonances and widths must consistently 

fit spectra for all angles. The high background and low 

counts obtained at 120° made that data comparably erratic and 

2 

the results less reliable despite the relatively low x • 
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Figures 1 through 4 show the experimental inelastic 

spectra with the fitted total background and the individual 

resonances superimposed. Figures 5 through 8 show the 

corresponding spectra with the radiation tail subtracted 

revealing the resonance structures more clearly. 

. 2 

Inelastic form factors squared, F , for each resonance 
are computed by the fit program, as reported in Ref. [(33)], 
using 

A . 

|F ± | 2 = (~~) |F el | 2 (IH-2) 

el 

where F^ = inelastic form factor for resonance 

A^ = area under resonance 

A e ^-= area under elastic peak 

F . = elastic form factor, 
el 

B(EL) values are calculated with the form factor squared as 
calculated with DWBA, where 

F. 2 

B (EL) = . (III-3) 

f dwba 

Tables II through V give the form factor squared and B value 
for each resonance at each angle. Figures 9 and 10 give 
the angular dependence of the squared form factor for the 
electric and the magnetic multipolarities with the curves 
normalized to the same maximum height. These curves were 
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used in making multipolarity assignments as shown in 
Figures 11 through 20. 

Table VI presents the compiled resonances determined, 
their energies, widths, B values and multipolarity assign- 
ments. Energy weighted sum rule per cent and single 
particle unit calculations are presented for each resonance 
in Table VII. 

D. ERROR ANALYSIS 

As line shapes for each of the resonances were fit to 

the four spectra of 75°, 90°, 105° and 120°, obvious 

variations in energy location and widths were noted. These 

variations were limited to approximately 200 KeV in both 

2 

energy and width while still maintaining reasonable x ' s * 

A definite correlation was noted between radiation tail 
subtraction changes in the fit program and the resulting 
transition strength form factor squared outputs. Table VI 
shows the largest deviation from the average values of 
resonance energies, widths and form factors of the four 
angles for each resonance. The reasons for these relatively 
large errors including statistical fluctuations, equipment 
functional variance and operator error. Warshawsky and 
Weber [Ref. (31) ] report the various LINAC equipment and 
operator contributions to errors. The repeatability of the 
90° and 120° spectra both indicate considerable improvement 
in the stability of the LINAC, as compared with Ref. 31. 

To determine the statistical error of the average reduced 
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transition matrix element, a weighted B (EL) -value was 
computed using each of the four spectrum values. 

Z B i (EL)/(AA i ) 2 

B (EL) = — 

Z 1/ (AA. ) 2 
i 1 

where: i = 75°, 90°, 105°, 120° 

AA^ = peak area statistical error from the fit 
program. 

The statistical error of this weighted average B (EL) -value 
is then 

B (EL) 

AB(EL) = . 

Z l/(AA i )" : 
i 

These errors were consistently less than 1% of the computed 
B (EL) values, much less than the variance of the B^ . Hence, 
this statistical error was mainly used to calculate B. 

The total error given in Table VI is a rounded combination 
of the statistical error and a best estimate obtained from 
fluctuations in the value of the form factor squared 
experienced during the fitting process . These errors ranged 
from 10 to 30%. 
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TABLE II* RESONANCES AT 75° SCATTERING ANGLE 



E 

(MEV) 


a 

-l 

(FM ) 


XL 


FCPM 

FACTOR 

SGUAPEC 


E VALUE 
(* ) 


6.60 


.553 


E2 


-4 

2.37X10 


2.11X10 






Ml 




2.55X10 


E • G9 


.549 


52 


-4 

2.76X1C 


2.62X10 






Ml 




4 .06X10 


5.£0 


.544 


E2 


-5 

8.77X1C 


1. 15X10 


11.15 


.540 


E2 


-5 

6.63X1C 


£.70X10 






EO 




1. 54X10 


12.60 


.536 


E3 


-5 

3.94X10 


6. 11X10 






M2 




2. 15X10 


12.50 


.534 


M2 


-5 

7.71X10 


6.27X10 






53 




1. 19X10 


14.60 


.526 


E 2 


-4 

2.53X1C 


2.22X10 


16.60 


.526 


El 


-4 

7.19X10 


1.57X10 


2 7.70 


.500 


E2 


-4 

4.71X1C 


6.26X10 













2 2L 

* UNITS: e FM FOR ELECTRIC 



2 2L-2 

e FM FOP MAGNETIC 



2 

0 

2 

0 

2 

I 

2 

| 

3 

1 

I 

A 

2 

I 

2 
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TAELE III. RESONANCES AT 90° SCATTERING ANGLE 













E 


Q 


XL 


FORM 


E VALUE 








FACTOR 






-1 




SCLAREC 


<* ) 


( ME V ) 


(FM ) 














-4 


2 


in 

r- 

• 

>4J 


.640 


E2 


1.60X1C 


2.52X10 










0 






Ml 




4.57X10 










2 


o 

c\j 

• 

»JU 


.635 


E2 


2.C5X1G 


3. 19X10 










0 






Ml 




5.62X10 








— c; 


1 


1C. 20 


.628 


E2 


4.10X1C 


6.21X10 








-5 


1 


11.25 


.625 


E2 


5.32X1C 


6.CCX1C 










2 






EO 




1 .44X10 






• 


-5 


3 


12.50 


.621 


E3 


5.33X1C 


5.63X10 










1 






M2 




3.53X10 








-4 


1 


12.60 


.617 


M2 


1.02X10 


6.67X10 










4 






E3 




1. 12X10 








-4 


2 


15.00 


.612 


E2 


3.04X10 


4.29X10 








-4 


1 


16.60 


.607 


El 


3.84X10 


1.97X10 








-A 


2 


2 6.10 


.572 


E 2 


3.42X10 


4.60X10 




_ 






ii - 



2 2L 

* LN ITS : e FM FOR ELECTRIC 



2 2L-2 

e FM FOR MAGNETIC 
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TABLE IV. RESONANCES AT 105° SCATTERING ANGLE 



E 

(MEV) 


Q 

-1 

(FM J 


XL 


FCPM 

FACTOR 

SCLAPEC 


E VALUE 
(* ) 


6.69 


.718 


E2 


-4 

1.16X10 


2.70X10 






Ml 




6.66X10 


S.C5 


.713 


E 2 


-4 

1.41X10 


2.20X10 






Ml 




1.C2X10 


1C. 15 


.705 


E2 


-5 

6.08X10 


1.21X10 


11.35 


.700 


E2 


— c 

4.24X10 


6.66X10 






EO 




1.60X10 


12.45 


.696 


E3 


-5 

3.57X10 


3. 57X10 






M2 




2.27X10 


12.60 


.691 


M2 


-5 

9.29X1C 


6.15X10 






E3 




9.29X10 


14.65 


.687 


E2 


-4 

2.37X10 


4.59X10 


16.60 


.680 


El 


-4 

1.71X10 


1.57X10 


27.80 


.639 


52 


-4 

3.06X10 


4. 61X10 













2 2L 

4 UNITS: e FM FOR ELECTRIC 



2 2L-2 

e FM FOR MAGNETIC 



2 

0 

2 

1 

2 

1 

2 

3 

1 

1 

2 

2 

1 

I 

2 
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TABLE V. RESONANCES AT 120° SCATTERING ANGLE 



E 

(MEV) 


Q 

-1 

(FM ) 


XL 


FORM 

FACTOR 

SCLAREc 


E VALUE 
(* ) 


6.73 


.783 


E2 


-4 

1.25X10 


2 

5.45X10 






Ml 




1 

1. 03X10 


8.00 


.778 


E2 


-4 

1.60X10 


2 

6.41X10 






Ml 




1 

1.42X10 


5. 5 0 


.770 


E 2 


-5 

3.10X10 


2 

1. 15X10 


11.10 


.765 


E 2 


_ c 

2.02X10 


1 

7. 18X10 






EO 




2 

1.31X10 


12.30 


.759 


E3 


-5 

5.08X10 


3 

5. 25X10 






M2 




1 

3 • 6 £ X 10 


13. £0 


.753 


M2 


_ G 

9.17X10 


1 

6 • SOX 10 






E3 




S. 55X10' 


IS. 00 


.748 


E2 


-5 

9.83X10 


2 

2.00X10 


16.60 


.741 


61 


-4 

1.08X10 


1 

1.57X10 


27.80 


.695 


E2 


-4 

2.71X10 


2 

5. 52X10 



2 2L 

* UMTS: e FM FOR ELECTRIC 



2 2 L-2 

e FM FOR MAGNETIC 
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TAdE VI. COLLECTIVE RESONANCES ANC I C ENT I F IC AT IONS 
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FCP MAGNETIC 



TABLE VII. SUN RULES ANC SINGLE PARTICLE TRANSITION 



STRENGTHS FOR 



E 

(M5V ) 


XL 


EWSR 

% 


6.69 


E2 


7.711. 


6.09 


E2 


11.0+2. 


1C. Cl 


E2 


4 • 3±0 • 


11.21 


E2 


3.7+1. 


12.46 


53 


3.5+0. 


1 -.62 


M2 




14.66 


E2 


24.6+4. 


16.60 


El 


100.6+0. 


27.68 


E2 


46.2+13. 



Y 



e 

(SPL ) 


E 

-1/3 

(A MEV) 


2.5 


30 


2.9 


36 


0.9 


45 


0.7 


50 


1.6 


56 


35.4 


61 


3.6 


66 


5.1 


74 


4.6 


125 



89 

6 

2 

4 

I 

7 

5 

1 

9 
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COUNTS PER JLXC (xIO 
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7 5° Spectrum with background 
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COUNTS PER UC (xIO 
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FIGURE 4. 



12 0° Spectrum with background 
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Form Factor Squared 




FIGURE 9. 
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Form Factor Squared 




FIGURE 10. 



Ml • M2 



Form Factors in comparison 



Form 




FIGURE 11 . Experimental inelastic form 

factors for state at 6.69 MeV 



Form 




FIGURE 12 • Experimental inelastic form 

for state at 8.09 MeV 
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FIGURE 13 . Experimental inelastic form 

factors for state at 10.01 MeV 
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FIGURE 14. 



Experimental inelastic form 
factors for state at 11.21 MeV 
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Form 




FIGURE 15 . Experimental inelastic form 

factors for state at 12.46 MeV 
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Form 




FIGURE 16 . Experimental inelastic form 

factors for state at 13.63 MeV 
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Form 




factors for state at 14.86 MeV 
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Form 




FIGURE 18. Experimental inelastic form 

factors for state at 16.60 MeV 
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Form 




FIGURE 19 . Experimental inelastic form 

factors for state at 27.85 M eV 
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IV. DISCUSSION 



A. COLLECTIVE RESULTS 

A total of nine transitions were determined to give 

consistent fits to the spectra at 75°, 90°, 105° and 120° 

2 

(See Table VI) . The x error per degree of freedom ranged 
from 0.86 at 105° to 1.08 at 75°. Observed variations in 
excitation energy and widths for the nine resonances were 
approximately 200 keV from their numerical averages. With 
few exceptions, the form factor values extracted from the 
fitting program compared favorably to DWBA calculations for 
the selected resonances. Some inconsistencies were noted in 
the extracted form factors for the 120° data. The relatively 
high background and possible transverse contributions to the 
cross section were contributing factors to these inconsistencies. 

The iterative process of fitting the spectra was started 
by fixing the excitation energy, width and ratio of elastic 
to inelastic peak height for the El giant dipole resonance 
[Refs. (28) and (29)] at E =16.6, F = 4.0 MeV and a peak 

X 

height which resulted in a reduced transition strength of 
2 2 

B(El)»l9.7e fm. Once this was done the remaining data were 
fitted by visually determining possible resonance excitation 
energies and widths and allowing the fit program to search 
for appropriate heights. Throughout the early stages of 
data analysis, repeated attempts were made to fit the 
previously reported E2 giant quadrupole resonance with 
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E x 13.8 and r 3.2 [Ref. (25)] to the data. Although 
2 

values of x <1.1 were obtained in these attempts, the 
data showed a very distinct structure at E =13.6 MeV 

A 

which could not be satisfactorily fitted in conjunction 

with the reported E2 values. The inclusion of the resonance 

at E = 13.6 MeV with V =1.2 MeV resulted in movement of 

the giant E2 resonance to E^ = 14.8 and T = 3.0 to correctly 

fit the data. The 90°, 105° and 120° spectra strongly 

supported the necessity for such a resonance and although 

its presence was less apparent in the 75° data, the 

resultant fit showed its existence. Once the fitting program 

was completed in the primary data sets, the resultant 

excitations were applied to the two alternate 90° data sets 

2 

and the one extra 120° data set. Excellent (x < 1.0) 
results were immediately obtained. 

B. EXCITATIONS AND CHARACTERISTICS 
1. 6.69 MeV 

This state is a very well defined narrow line of 
T = 0.97 MeV which was given the assignment of an E2 bound 
state. Its width was twice the resolution width. The line 
structure was quite smooth in the 75°, 90°, and 105° dhta 
but became more erratic in the 120° data. The latter 
behavior can be explained by the low count rate and high 
background and the increasing transverse contributions to 
the scattering cross section at the backward angle. Using 
this reduced transition strength, the resonance accounts 
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for 7.7 ± 1.6% of the EWSR and 2.5 Weisskopf units for the 
E2 assignment. Visual inspection of form factor plot, 

Figure 11, favors the E2 assignment. 

2. 8.09 MeV 

This resonance is the second well defined line in 
the bound state region. The state has a width of 
T = 1.22 MeV and was also given an E2 assignment. Again 
the line structure in the data was quite smooth and easy 
to fit at 75°, 90°, and 105°. The assignment was made 
based on these angles and the resultant B(E2) value accounts 
for 11.0 ± 2.2% of ESWR and 2.9 Weisskopf units. 

3. 10.01 MeV 

-1/3 

The assignment of E2 at 45A was based on 75°, 

105° and 120° data. Considerable difficulty was encountered 
in attempting to fit the energy range from 10 MeV to 12 MeV 
due to the relatively slow count rate in comparison to the 
twin peaks at 6.69 MeV and 8.09 MeV and the giant resonance 
region. The resonance at 10.01 MeV was consistently erratic 
in all three 90° spectra since the number of counts was very 
low at that energy. The assignment of E2 accounts for 
4.3 ± 0.4% of the EWSR and 0.9 Weisskopf units. 

4. 11.21 MeV 

The analysis of the state at 50A MeV favored an 

E2 assignment. Since the DWBA form factors for E0 and E2 are 
very similar in the momentum transfer region of this study, 
both assignments were considered. The E2 assignment accounts 
for 3.7 ± 1.1% of the EWSR and has a strength corresponding 
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to 0.7 Weisskopf units. However, an EO assignment exhausts 
5.8% of the available monopole strength in Yttrium. 

5. 12.46 MeV 

The multipolarity assignment of the resonance at 

-1/3 

56A ' remains ambiguous. The extracted form factor values 

agreed most closely with an M2 assignment although an E3 
assignment served almost as well. The M2 assignment 
displayed a strength of 17.5 Weisskopf units. An assignment 
of E3 exhausted 3.5 ± 0.7% of the EWSR and yielded a strength 
of 1.6 Weisskopf units. 

6. 13.63 MeV 

-1/3 

The existence of this resonance (at 61A ) was 

strongly supported by all of the data. It was given an 

assignment of M2. The relatively narrow structure appeared 

in each spectrum, including the three separate 90° runs and 

the two separate 120° runs. Its presence was less obvious 

in the 75° spectrum but inclusion of the resonance improved 

8 9 

the fit. The E2 giant resonance in Y has been reported 

at 13.8 ± .2 MeV with a width of 3.2 from (p,p') work 

90 

[Ref. (25) ] . Electron scattering experiments in Zr have 
reported the E2 GDR at E =14.0 MeV and V = 4.8 MeV 
[Ref. (4)]. Therefore, the assignment of the relatively 
narrow line at this energy proceeded cautiously. 

The extracted form factors compared favorably to 
both a M2 and an E3 assignment. An E3 assignment would have 
exhausted 7.7 ± 1.5% of the EWSR with a strength of 3.2 
Weisskopf units. The M2 assignment had a reduced transition 
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strength correlating to 35.4 Weisskopf units. The electron 
scattering experiments should have shown the contribution of 
such a strong E3 mode if it were indeed present. However, 
in the latter experiment, an M2 resonance would have been 
obscured by the E2 because M2 and E2 display a similar 
dependence on q at very forward angles where the transverse 
character of a magnetic transition does not matter very much. 
7. 14.86 MeV 

-1/3 

The assignment of an E2 giant resonance at 66A 

8 9 

compares to some extent with the (p,p') experiments in Y 

reported in Ref. [(25)] and with the (e,e') experiments in 
90 

Zr reported in Ref. [(4)]. The shift in excitation energy 

is a direct result of fitting the required resonance at 

-1/3 

13.63 MeV. The data readily support the 66A energy 

position and width of 3.0 MeV. Extracted form factors 
compared well with DWBA calculations for an E2 resonance 
and the assignment exhausts 24.6 ± 4.9% of the EWSR with a 
strength of 3.6 Weisskopf units. This relatively low 
exhaustion of EWSR compared very well with the results of 
the (p,p') experiments but was less than half the corre- 
sponding value reported for the (e,e') and the (a, a') 

9 0 

experiments in Zr. The latter observation prompted 
extensive attempts to fit the data with an E2 cross section 
that would account for about 50% exhaustion. The surprising 
result was that 90°, 105° and 120° data would support such 
a cross section with E x = 14.0 MeV and f = 4.8 MeV and an 
exhaustion of 53% E2 EWSR. However, the 75° data, which 
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displayed a much smoother character in this area, would not 

accept such a large cross section. The largest E2 cross 

section able to fit the 75° data exhausted 35% EWSR and 

2 

resulted in a poor fit (x =1.4). Repeated attempts to 
match the E2 cross sections through the four angles led 
irrevocably to the results reported in this paper. For a 
comparison of these different fitting attempts see Figures 6 
and 20. It was also noted that (a, a*) experiments reported 
in Ref. [(24)] describe the E2 giant resonance at E = 14.5 

X 

MeV, T = 4.0 MeV with an E2 EWSR exhaustion of 54%. Again 
no resonance of this description could be coaxed to properly 
fit all the available spectra. 

8. 16.6 MeV 

-1/3 

The resonance at 74A is the known El giant dipole 

resonance reported from photonuclear work [Refs. (28) and (29)] 

and compares well with the reported position, width and 

90 

strength of the El GDR in Zr electron scattering experiment 

[Ref. (4)]. Considering the paucity of experimental data in 
89 

Y, this resonance was used exclusively as a known, fixed 

feature for the data fitting process. The results obtained 

on all other fitted resonances are therefore greatly dependent 

upon the correct assignment of E = 16.6, T = 4.1 and reduced 

2 2 

transition strength of 19.7 (e MeV fm ). These values were 
observed to fit the data very well and resulted in exhausting 
100.6% of the EWSR with a strength of 5.1 Weisskopf units. 

The excellent results obtained from fixing this resonance 
lent confidence to the assignments of the remaining resonances. 
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Some evidence exists in the 75° data for the excitation of 



the T > [Ref. (34)] part of the El giant resonance at 
E x = 20.3 MeV, but attempts to fit T > were unsatisfactory 
and the remaining angles offered no support for this 
resonance . 

9. 27.92 MeV 

-1/3 

The resonance at 12 5A was given an E2 assignment 

and compares well with the predicted E2 isovector resonance 
of 135A [Ref. (18)]. This assignment exhausted 

46.3 ± 13.9% of the isovector E2 EWSR and the transition 
strength accounted for 4.6 Weisskopf units. 
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V. CONCLUSIONS 



The spectrum of inelastically scattered electrons from 
89 

Y was studied in an excitation energy range from 6.1 MeV 
to 38.0 MeV. Five previously unreported states were 
observed below the (y,o) threshold. Four giant resonances 
were observed at 13.63, 14.86, 16.6 and 27.92 MeV. Of 
these, the 13.63 MeV M2 and the 27.9 MeV E2 isovector 
resonances have not been previously reported. Multipolarity 
assignments were made as discussed in Section IV. B. 

The 16.6 MeV resonance was identified as the widely 
studied El giant dipole resonance and was used as a starting 
value for data analysis. Although the assignments of the 
M2 at 13.63 MeV and the E2 giant quadrupole at 14.8 MeV 
do not support existing studies in this and neighboring 
elements [Refs. (25), (4) and (24)], the assignments were 

considered positive. Previous (e,e') studies of a N = 50 
nucleus were conducted at forward scattering angles 
[Ref. (4) ] and therefore would not be expected to distinguish 
the M2 from the E2 cross section. Transverse magnetic 
contributions may also have been present in the 6.69 and 
8.09 MeV states at 120°. Further study in the backward 
angles should be made with better statistics. 

The E2 resonance at 27.92 MeV was observed as predicted 
in Ref. [(18)] and is believed to be the isovector E2 giant 
resonance . 
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